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NOTATION  (CONT.) 
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CG 
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3C  /3 (rb/2  V  ) ,  per  radian 
n  e 

2 

weight  coefficient,  W/J  PV  S 

2 

side  force  coefficient,  Y/$  PV  S 
SC^/SE,  per  radian 

3C  /3(pb/2  V  ),  per  radian 
y  e 

3C_/3(rb/2  V  ),  per  radian 
T  e 

mean  aerodynamic  chord,  1.50  m  (4.91  ft) 
centre  of  gravity 
drag  force,  N 

2 

acceleration  due  to  gravity,  m/sec 

distance  of  C.G.  from  wing  leading  edge,  expressed  as  a 
fraction  of  c 

altitude,  m  (2000  ft  in  this  study) 

2 

moment  of  inertia  about  X  body  axis,  1395  kg-m 
(1025.6  slug-ft2) 

2 

moment  of  inertia  about  Y  body  axis,  1480  kg-m 
(1088.1  slug-ft2) 

2 

moment  of  inertia  about  Z  body  axis,  2563  kg-m 
(1884.3  slug-ft2) 

2  2 
product  of  inertia,  123  kg-m  (90.4  slug-ft  ) 

distance  from  horizontal  tail  aerodynamic  centre  to  the 
aircraft  C.G.,  4.51  m  (14.8  ft) 

pitching  moment,  N-m 

aircraft  mass,  1318.2  kg  (90.06  slugs) 

yawing  moment,  N-m 

normal  acceleration,  g  units 

roll  rate,  rad/sec 

dynamic  pressure 

lateral  gust  gradient,  -(l/v  )  (dw  /dy) 

e  g 

pitch  rate ,  rad/sec 
yaw  rate,  rad/sec 

2  2 

wing  area,  16.26  m  (175  ft  ) 

3  2 

horizontal  tail  area,  3.59  m  (38.65  ft  ) 
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NOTATION  (CONT.) 


thrust,  N 

2 

thrust  coefficient,  T/J  PV  S 
airspeed,  m/sec 

horizontal  tail  volume  ratio,  S^/Sc  (0.67) 

aircraft  weight 

side  force,  N 

angle  of  attack,  radians 

angle  of  attack  due  to  w  gust,  (w  /V  ) 

g  e 

angle  between  thrust  vector  and  fuselage  ref.  line 
(-3.5  deg.) 

angle  of  sideslip,  rad 

angle  of  sideslip  due  to  v  gust,  (v  /V  ) 

g  e 

elevator  deflection,  positive  when  trailing  edge  is  down, 
deg. 

damping  ratio 
roll  angle,  deg 
climb  angle,  deg 

3  2 

air  density,  kg/m  (0.0022  slug/ft  at  2000  ft) 
pitch  angle,  deg 

input  gust  intensity,  m/sec  (10  ft/sec  in  this  study) 

standard  deviation  of  aircraft  velocity  response,  ft/sec 
(fps  units  in  program) 

standard  deviation  of  angle  of  attack  response,  deg. 

standard  deviation  of  pitch  rate  response,  deg/sec 

standard  deviation  of  pitch  angle  response,  deg 

standard  deviation  of  normal  acceleration  response, 
fractional  g 

standard  deviation  of  aircraft  sideslip  response,  deg 
standard  deviation  of  aircraft  roll  rate  response,  deg /sec 
standard  deviation  of  aircraft  yaw  rate  response,  deg/sec 
standard  deviation  of  aircraft  roll  angle  response,  deg 
standard  deviation  of  aircraft  yaw  angle  response,  deg 
time  constant,  seconds 
yaw  angle,  deg 
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Subscripts 

e  Parameter  or  coefficient  at  triamed  condition 

g  Refers  to  gust 

A  dot  over  a  symbol  signifies  a  derivative  with  respect 
to  time 
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1.  INTRODUCTION 

The  work  reported  here  is  part  of  a  task  to  determine 
the  response  spectra  of  various  aircraft  due  to  gust  turbulence.  It 
was  done  for  the  Advanced  Engineering  Laboratory,  Salisbury,  which 
requires  data  for  use  in  the  design  of  inertially  stabilised  airborne 
control  systems.  This  work  has  also  served  as  an  introduction  to  the 
concepts  of  aircraft  stability  and  control,  to  augment  the  author's 
instrumentation  and  control  background. 

The  aircraft  configuration  which  is  the  subject  of 
this  Memo  is  broadly  representative  of  high  wing,  general  aviation 
monoplanes.  Specifically,  the  aerodynamic  data  used  in  the  mathematical 
models  have  been  obtained  from  Reference  2  and  Reference  3.  Estimates 
for  centre  of  gravity,  weight  and  inertias  have  been  made  using 
information  in  Reference  3  and  Reference  4  as  a  guide. 

This  Memo  outlines  the  procedure  used  to  calculate  the 
various  stability  derivatives  and  tabulates  all  the  input  data  to  the 
GUSTR  program  described  in  Reference  1 .  Output  frequency  responses 
and  spectra  results  are  included  and  a  table  of  calculated  responses 
to  an  input  gust  intensity  of  10  ft/sec  is  given. 

2.  LONGITUDINAL  PARAMETERS 
2.1  Trim  Determination 


The  parameters  for  three  CG  positions  of  15%,  25%  and 
35%  of  mean  aerodynamic  chord,  c,  were  evaluated.  These  three  positions 
were  chosen  on  the  following  basis. 

The  moment  reference  for  the  model  in  Ref.  2  was  0.2  m 
from  the  wing  leading  edge,  or  13.34%  c.  Ref.  2  also  indicated  that  the 
model  was  neutrally  stable  with  the  CG  as  far  aft  as  40%  c.  Therefore, 
three  CG  positions  of  15%,  25%  and  35%  c  were  chosen  for  this  study, 
together  with  two  speed  cases  of  120  knots  and  160  knots.  This  resulted 
in  a  total  of  six  sets  of  trim  angle  of  attack,  a,  elevator  angle,  6  , 
and  thrust  coefficient,  T",  together  with  six  corresponding  sets  of 
longitudinal  stability  derivatives.  An  altitude  of  2000  ft  and  zero 
climb  angle  have  been  assumed  for  all  cases. 

The  trim  values  of  a,  6  ,  and  T*  were  obtained  by  first 
tabulating  the  C',  C',  C'  curves  in  fig.  7(a)  ,C(b),  (c)  of  Ref.  2  for 

IP  L  D 

06  <  a  <  10° 

-15*  <  6  <  5° 

e 

0°  <  T'  < 
c 


0.3 
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This  file  of  540  data  points  was  used  in  a  least  squares  fitting 
program  (Fig.  1)  to  provide  coefficients  for  C^«  and  in  the 
following  (not  necessarily  linear)  representation 
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Note  that  the  representations  of  forces  and  moments,  C^,  and 
C^,  given  in  Equations  (1)  to  (3)  include  both  aerodynamic  and 
thrust  contributions.  Consequently  the  trim  equations  are  (see 
Fig.  2) 


CL-CWSin\=°  ••'(4) 

C'  +  C  cos  y  =  0  ... (5) 

D  W  e 

C'  =0  ...  (6) 


where  Cw  =  W/iPV  S,  Y  is  the  steady  climb  angle  and  C£,  and  C'  are 
given  as  functions  of  a,  6  and  T'  by  Equations  (1)  to  (3)  above. 

For  given  values  of  density,  p,  aircraft  weight,  W,  trim  speed,  V  and 
climb  angle,  Y  ,  Equations  (4)  to  (6)  constitute  a  set  of  non-linear 
algebraic  equations  for  trim  values  of  a,  6  and  T'.  The  solution  is 
obtained  in  program  TRIM  (Fig.  1)  and  the  risults  for  the  six  cases 
considered  are  given  ir.  Table  1. 

When  it  is  necessary  to  separate  out  the  aerodynamic 
component  of,  say,  from  the  thrust  component,  this  can  be  done  by 
noting  that  Equation  (2)  can  be  rewritten  (see  Fig.  2) 

a  a 

C  -T '  cos(a+a„)  =  b-  +  b, -a  +  b„.a  +  b,.c  +  b..6  +  b_.a£  +  b, .  T' 

Dc  T  01  2  3  e  4  e  5e6c 

or,  for  small  (a  +  a_) 

T 

C  *=  b-  +  b  a  +  b.a  +  b,6  +  b.£  +  bca£  +  (b,  +  1)  T' 

D  0  1  2  3e4e  5e  6  c 


... (7) 
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where  is  now  the  aerodynamic  component  only.  Similarly,  in 

Equation  (1)  c£  can  be  replaced  by  C  +  T'  sin(a+a  )  which  is 

approximately  equal  to  C  (Aerodynamic)  for  small  (a+a  ) . 

L  T 

The  coefficients  for  Equations  (1)  to  (3)  obtained 
from  the  least  squares  fitting  program  are  listed  in  Table  2.  A 
sample  check  on  their  accuracy  was  made  by  verifying  that  Equations 
(4)  to  (6)  were  satisfied  when  the  trim  values  of  a,  6  and  T'  were 

Q  C 

substituted  in  these  best  regression  solutions. 

2 . 2  Aircraft  Weight  and  Inertias 

The  CG  envelope  in  Reference  4  gives  a  ratio  of 

Airplane  Weight  (Most  forward  CG  position)  =  0.81 
Max.  Take  off  weight 

Using  the  same  factor  on  the  aircraft  of  reference  2 , 
with  maximum  take-off  weight  =  3600  lbf  gives 

W(fwd  CG)  =  2916  lbf. 

Therefore  a  figure  of  2900  lbf  for  aircraft  weight 
was  used  throughout  this  study,  being  the  smallest  maximum  allowable 
weight  within  the  CG  envelope  considered. 

The  inertias  Ix,  Iy,  Iz  and  Ixz  were  taken  from 
Reference  3  and  were  not  corrected  for  the  different  CG  positions 
of  this  study.  These  values,  together  with  those  of  the  pertinent 
geometric  characteristics  used  in  this  study,  are  listed  in  the 
notation. 


2 . 3  Determination  of  longitudinal  Stability  Derivatives 

The  combination  of  wind  tunnel  results  (Ref.  2)  and 
flight  measurements  (Refs.  3  and  4)  was  the  source  of  some  derivatives, 
while  others  not  directly  measured,  were  derived  using  the  formulae 
of  Reference  5.  Note  that  the  static  derivatives,  e.g.  CL  ,  C^  (in 
the  lateral  case)  etc.  were  obtainable  from  both  wind  a  6 
tunnel  and  flight  tests.  Dynamic  derivatives,  e.g.  /fc  +  C^  N,C£ 

(in  the  lateral  case)  etc.  were  only  obtainable  fromV  mq  md /  lp 

the  flight  tests.  Also,  C  and  C  could  not  be  separated  from 

m  m. 

q  a 

each  other  in  the  flight  test  results,  therefore  this  was  an  instance 

where  C  was  first  calculated  and  C  then  followed  from  the 

m  m. 

q  a 

subtraction  of  C  from  the  combined  flight  test  result  of 

q 


( C  +  C  )  . 
V  m  m.  ) 


2.3.1  Static  derivatives 


Parameters  CT  ,  C  ,  C  ,  C  ,  C  ,  C,  ,  C  ,  C  , 

L  D  L  D  m  L.  m,  D . 
eeaaao  o  o 

e  e  e 

C  ,  C  were  evaluated  directly  from  the  least  squares  representations 

°T '  mT - 
c  c 

of  Eqns.  (1),  (3)  and  (7).  E.g.  From  Eqn.  (7) 

3C_ 


"D  3a 

at 


b  +  2b.. a,  .  - 
1  2  Trim 


2.3.2  The  V  derivatives 

The  aircraft  has  a  constant  speed  propeller, 
therefore  (Ref.  5) 


-3C_ 


The  derivative  CL  has  an  expression  typical  of  the  V 
derivatives,  where  (Ref.  5)  v 
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_ L 
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T  ST" 
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fl 

3M 


f  0  below  0.6M  (negligible  compressibility  effects) 


7  0  (negligible  aeroelastic  effects) 


and 


ffl 

3T ' 


-  0  (from  the  Least  Squares  representation  of  CL) 


Thus  C. 


0. 
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Similarly,  the  first  two  terras  disappear  in  the  expression  for  Cc 
and  ,  leaving 


c  .  c 

T  D_ 
v  T' 


C  .  C 
t  m 
v  T " 


2.3.3  The  q  derivatives 

C  =  -  2C  (h  -  h  >  +  (2a  V  )  (Ref.  5) 

Li  Li  O  t  H 

q  a 

where  h  =0.75  for  subsonic  flow 

o 

VH  it  st  =  0.67 
Sc 

=  14.8  ft  (scaled  from  fig.  2,  Ref.  2) 

Sfc  =  30.65  ft*  (Ref.  7) 

=  3.72/rad  (Calculated  from  Ref.  8,  knowing  the 
tail  aspect  ratio) 

Therefore  for  Case  1,  with  CG  at  15%  c 
h  =  0.15. 

Thus  C  =6.40+4.99 

L< 

q 


C_  =  11.39/rad. 

Ij 

q 


Also 


cra  *  £m  -  2CL  <h  '  -  2at  «t 

q  q  a  - 


(Ref.  5) 


For  subsonic  flow,  C  =0 


h  =  0.5 


Therefore  fore  Case  1 


C  =  0  -  1.31  -  15.04 

m 

q 

C  =  -  16.35/rad. 

rr. 

q 


2.3.4  The  d  derivatives 


Ref.  3  gives  C  +  C  -18.5  at  C,  =0.36 
v  m  m.  L 

q  a 


Hence  C  =  -2.15/rad  for  Case  1. 

m. 

a 


and  C.  =  -  C  .  C  =  0.713/rad/sec  for  Case  1. 

L.  m. 

a  a 


A  summary  of  the  Longitudinal  derivatives  is  given  in 

Table  3. 


3.  LATERAL  DERIVATIVES 


The  derivatives  C 


n 


y« 


P  P  P  r  r 
were  all  derived  from  the  flight  test  results  of  Reference  3.  The 
'maximum  likelihood1  results  were  chosen  and  interpolated  values  of 
the  derivatives 
were  obtained. 

wind  tunnel  results  of  Reference  2,  were  used  as  a  cross-check.  C 
agreed  within  1%.  c  ,  after  suitable  correction  for  CG  position, 


for  C  =0.36,  corresponding  to  120  knots  at  2000  ft. 
The  6  derivatives,  which  were  also  available  from  the 


6 


was  found  to  be  10%  higher  and  C»  was  found  to  be  7.5%  higher  in  the 

e 

wind  tunnel  results.  A  listing  of  the  lateral  derivatives  used  is 
given  in  Table  4. 


The  lateral  gust  response  was  computed  for  V  =  120  knots 
at  2000  feet  with  CG  located  at  0.15  c.  Other  cases  can  be  run  as 
required. 


given  at  a  particular 


The  method  of  correcting  , 

CG  location  (e.g.  0.28  c  in  Ref.  3),  for  a  different  CG  location, 
(e.g.  0.15  c  in  the  lateral  case  computed),  is  as  follows 


C  =  C  -  (0.26  -  0.15) .c  .  C 

n .  n, 

C (0 . 15  c)  c(0. 28  c) 


b 


E.g.  for  C  =  0.050/rad  and  C  =  -0.59/rad 

n,  y  „ 

S (0 . 28  c)  b 


C  =  0.060/rad. 

nB 

<0.15  c) 


4 .  RESULTS 


The  results  of  all  longitudinal  and  lateral  cases  run 
are  tabulated  in  Tables  5  and  6.  The  computed  frequency  responses 
and  output  spectra  of  the  GUSTR  program  (Ref.  1)  are  shown  in  Figs. 

3  to  22.  The  von  Karman  gust  spectrum  was  used,  with  an  o  input 
to  the  longitudinal  model  and  both  6  and  pg  inputs  to  theglateral 
model  (Ref.  1).  All  output  variancegcalculations  in  Table  5  were 
based  on  an  input  gust  intensity  of  o  =10  ft/sec.  and  the  formula 
used  was  (Ref.  1)  * 

2 

o 

— —  =  4.6  X  (area  under  curve)  X  (scale  factor) 

a 

g 


Note  that  the  computer  program  uses  radian  measure  and  f.p.s.  units 
in  the  output  plots. 

Areas  under  the  output  spectra  curves  were  measured 
with  a  planimeter  in  the  range 


4  4  Log  Cj  {  - 


1 
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corresponding  to  a  frequency  range  of 


0.0032  Hz  -*  3.2  Hz  for  the  120  knot  cases 

and  0.0043  Hz  -*•  4.3  Hz  for  the  160  knot  cases. 

The  theoretical  range  of  validity  is  at  '  =  0 . 2 

c 

i.e.  at  Log  ~L  =  -1.39. 

However,  for  comparing  the  relative  trends  between  the  various  cases, 
extending  the  cut-off  limit  to  Log  fi  =-1.00  (as  done  in  this  exercise) 
should  not  introduce  significant  error. 

5.  DISCUSSION 


The  results  are  generally  as  expected,  especially  the 
variation  of  natural  frequencies  and  damping  ratios  of  the  different 
modes  with  speed,  V  ,  and  CG  location.  It  is  worth  noting  that  at 
both  rearward  CG  locations  of  0.35  c  (Cases  3  and  6),  the  Short 
Period  Mode  has  real  roots. 


The  output  spectra  most  likely  to  be  of  relevance  to 

the  designer  of  inertial ly  stabilized  airborne  control  systems  are 

the  attitude  and  rate  spectra,  viz.  6  and  q  spectra  in  the  longitudinal 

cases  and  V,  p  and  r  in  the  lateral  cases.  As  shown  in  Table  5, 

the  deviation  of  pitch  rate  response,  o  ,  decreases  strongly  as  the 

CG  shifts  rearwards  while  also  decreases  somewhat,  particularly 

at  160  knots.  There  is,  at  the  same  time,  only  a  slight  increase  in 

variance  of  the  normal  acceleration  response ,  o  .  Therefore  from 

n 

z 

the  inertial  platform  designers  standpoint,  a  CG  location  of  0.35  c 
would  be  desirable.  However,  a  more  practical  figure  of  CG  at  0.30  c 
may  be  a  better  compromise  between  acceptable  handling  qualities  and 
inertial  platform  design  requirements. 


Of  the  lateral  responses  in  Table  6,  the  dominant 
contribution  of  p  to  the  roll  response,  o  ,  is  due  to  the  large 
value  of  the  rollgdamping  derivative,  C^  .  For  all  other  response, 


the 


S  gust  is  dominant, 
g 


All  the  results  presented  here  are  for  the  aircraft 
’stick  fixed'  response.  Reference  4  suggests  that  the  'stick  free' 
response  of  these  aircraft  would  not  be  significantly  different  due 
to  control  linkage  stiction  and  friction. 
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6.  CONCLUSION 


This  Memo  has  outlined  the  steps  taken  to  obtain  the 
longitudinal  and  lateral  rigid-aircraft  gust  response  of  a  light 
single-engmed  high-wing  aircraft  configuration  for  a  few  specific 
cases  of  speed,  CG  location,  aircraft  weight  and  altitude.  Estimates 
of  all  required  static  and  dynamic  derivatives  are  provided  so  the 
results  for  any  other  combination  of  flight  conditions  or  CG  location 
can  be  readily  obtained. 
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TABLE  1  -  TRIM  CONDITIONS 


In  all  cases,  altitude,  h  =  2000  ft;  climb  angle,  Y  =  0°. 


Identifier 

,  Ve 
(Knots) 

CG  Position 
(%  c) 

3  (trim) 
(deg . ) 

5 

e (trim) 
(deg. ) 

Case  1 

120 

15 

1.87 

0.90 

Case  2 

120 

25 

1.70 

1.98 

Case  3 

120 

35 

1.56 

2.96 

Case  4 

160 

15 

-0.09 

2.67 

Case  5 

160 

25 

-0.11 

2.83 

Case  6 

160 

35 

-0.15 

3.07 

0.060 

0.060 

0.060 

0.060 

0.060 

0.060 


TABLE  2  -  BEST  REGRESSION  SOLUTIONS  FOR  C,  ,  C,  C 

L  D  in 


TABLE  3  -  LONGITUDINAL  STABILITY  DERIVATIVES 


CASE  1,  4 

CASE  2,  5 

CASE  3,  6 

CT  * 

D  D_ 

e  T' 

c 

0.057 

0.057 

0.057 

CL  (/rad) 
a 

5.330 

5.330 

5.330 

CD  (/rad) 
a 

0.118 

0.107 

0.098 

C  (/rad) 

in 
a 

-  2.029 

-  1.110 

-  0.292 

CT  (/rad) 

L6 

e 

0.797 

0.797 

0.797 

Cm  (/rad) 

6 

e 

-  2.239 

-  2.160 

-  2.080 

C  (/rad) 

U6 

e 

0.059 

0.059 

0.059 

T" 

c 

0.097 

0.099 

0.101 

CT 

V 

-  0.162 

-  0.162 

-  0.162 

CL 

V 

0 

0 

0 

CD 

V 

-  0.009 

-  0.009 

-  0.009 

Cm 

V 

-  0.016 

-  0.016 

-  0.016 

C  (/rad/sec) 

L 

q 

11.390 

10.320 

9.250 

Cm  (/rad/sec) 

q 

-16.350 

-15.710 

-15.280 

Cm  (/rad/sec) 
a 

-  2.150 

-  2.790 

-  3.220 

C  (/rad/sec) 

Li  • 

a 

0.713 

0.926 

1.068 

C,  =  C,,  *  0.36  for  cases  1,  2,  3 

L  W 


0.20 


for  cases  4,  5,  6 


FIG.  I  PROCEDURE  USED  FOR  EVALUATING  TRIM  CONDITION'S 


Power  Spectrum 


0 . 3  Hz  3  Hz 


LOG 


10 


(••) 


FIG.  5  PITCH  RATE  RESPONSE  TO  VERTICAL  GUST,  wQ ,  FOR 
V  =  120  KNOTS  y 


FIG.  7  NORMAL  ACCELERATION  RESPONSE  TO  VERTICAL  GUST,  wq, 
FOR  V  =  120  KNOTS 


FIG.  9  ANGLE  OF  ATTACK  RESPONSE  TO  VERTICAL  GUST,  wg , 
FOR  V  =  160  KNOTS 


B.  08 


FIG.  10  PITCH  RATE  RESPONSE  TO  VERTICAL  GUST,  w  FOR 
V  =  160  KNOTS  g 


SIDESLIP  ANGLE  RESPONSE  TO  LATERAL  GUST 
FOR  V  =  120  KNOTS 


15  YAW  RATE  RESPONSE  TO  LATERAL  GUST,  v  .  FOR  V 


10.00 


FIG.  18  SIDESLIP  ANGLE  RESPONSE  TO  LATERAL  GUST  GRADIENT,  P( 
FOR  V  =  120  KNOTS 


19  ROLL  RATE  RESPONSE  TO  LATERAL  GUST  GRADIENT 
V  =  120  KNOTS 


FIG.  20  YAW  RATE  RESPONSE  TO  LATERAL  GUST  GRADIENT,  p  , 
FOR  V  =  120  KNOTS 
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